Spin-transfer in difFusive ferromagnet- normal metal systems with spin-flip scattering 
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The spin-transfer in biased disordered ferromagnet (F) - normal metal (N) systems is calculated 
by the diffusion equation. For F1-N2-F2 and N1-F1-N2-F2-N3 spin valves, the effect of spin-flip 
processes in the normal metal and ferromagnet parts are obtained analytically. Spin-flip in the 
center metal N2 reduces the spin-transfer, whereas spin-flip in the outer normal metals Nl and N3 
can increase it by effectively enhancing the spin polarization of the device. 

PACS numbers: 72.25.Ba,75.70.Pa,75.60.Jk,72.25.Rb 
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A spin-polarized electric current flowing thmugh magnetic multilayers with canted magnetizations produces torques 
on the magnetic moments of the ferromagnetsEa The effect is inverse to the giant magnetoresistance, in which a current 
is affected by the relative orientation of the magnetization directions. The spin-current-induced magnetization torque 
arises from an interaction between conduction electron spins and the magnetic order parameter, transferring angular 
momentum between ferromagnetic layers, hence the name "spin transfer". The observed asymmetry of the switching 
with respect to the direction of current flow in the magnetization switching in cobalt layergjo™ is strong evidence 
that the spin-transfer dominates charge current-induced Oersted magnetic fields in mesoscopic small structures. Spin- 
transfer devices are promising for applications by the ability to excite and probe the dynamics of magnetic moments 
at small length scales. Reversing magnetizations with little power consumption can be utilized in current-controlled 
magnetk^Bjteajuap, .elgn,3|eijLt s . As a result the spin-transfer effect has already been the subject of several theoretical 
studies.fflil3eBlBye 

The torque can be formulated by scattering theory in terms of the spin-dependence of the refiection coefficients of 
the interface and the incoherence of spin-up and spin-down states inside the ferromagnet. This leads to destructive 
interference of the component of the spin-current perpendicular to the magnetization CLy|p_tiifj_&p:pmagnetic decoher- 
ence length, which is smaller than the mean free path for not too weak ferromagnets.l3allj'yOll3M In this paper we 
solve the spin-dependent diffusion equation for a multilayer system consisting of two reservoirs, three normal metal 
layers and two ferromagnetic layers, see Fig. |l|, generalizing the approach of Valet and FertEj to non collinear systems. 

We focus on relatively large systems in which the bulk resistance dominates. Interfaces play an essential role in 
transferring the torque, but are assumed not significantly to increase the total resistance. A typical interface resistance 
for e.g. Co/Cu is A-Rsurface ^ i^To^ . The corresponding typical bulk resistance for clean and dirty Co/Cu layer varies 
between Ai?buik = 0.01L[nm] iVtw? and 0.1i[nm] iVLw? (see for example Ref. |l^and Ref. |l^) where L is the length of 
the layers expressed in nm. Our approach is therefore quantitatively valid when the layers are thicker than 100 nm for 
pure samples and 10 nm for alloys which is reasonable for experimental fabrication, and furthermore reveal qualitative 
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FIG. 1: N1-F1-N2-F2-N3 multilayer system with non-collinear magnetizations. 
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effects of spin-flip relaxation process on the spin-torques for thinner layers. Related calculations of the torque and the 
magnetoresistance for submicron Co/Cu multilayers using the Boltzmann equation have been presented in Ref. |l^. 
Here we obtain analytical results based on the diffusion equation that reveal the main physical effects of spin-flip 
scattering in different parts of the multilayer on the spin transfer: Spin-flip scattering in the middle normal metal 
N2 reduces the spin-transfer whereas spin-flip scattering in the outer normal metals Nl and N3 can enhance the 
spin-transfer. 

The paper is organized as follows: in Section II we explain the averaging mechanisms of spin-transfer and the 
boundary conditions for the diffusion equation. The latter are formulated for a N1-F1-N2-F2-N3 multilayer system 
and solved analytically in the presence of spin-flip processes in the bulk layers in Sedjions III and IV. In the last 
section we summarize our conclusion. In Appendix A, magnetoelectronic circuit theory@li3 is shown to be consistent 
with the results from the diffusion equation in the absence of spin-flip scattering. 



II. DIFFUSIVE APPROACH TO MULTILAYER SYSTEMS 



Electron states with spins that are not collinear to the magnetization direction are not eigenstates of a ferromagnet, 
but precess around the magnetization vector. In three dimensions, a non-collinear spin current is composed by many 
states with different Larmor frequencies which average out quickly in a ferromagnet as a function of penetration depth. 
The efficient relaxation of the non-diagonal terms in the spin-density mairiy is equivalent to the suppression of spin- 
accumulation non-collinear to the magnetization in the ferromagnet BIlj'oEil This spin-dephasing mechanism does 
not exist in normal metals, in which the spin wave functions remain coherent on the length scale of the spin-diffusion 
length, which can be of the order of microns. In ballistic systems, the spin-transfer occurs over the ferromagnetic 

decoherence length \c — 1/ k\ — k\;, . In conventional ferromagnets the exchange energy is of the same order of 

magnitude as the Fermi energy and Ac is of the order of the lattice constant. The strongly localized regime in which 
the mean free path is smaller than the inverse Fermi wave- vector, (. < 1/kp, is not relevant for elemental metals. In 
conventional metallic ferromagnets i ^ I/^f, and the length scale of the spin-transfer Ac is necessarily smaller than 
the mean free path i and therefore is not affected by disorder. This argument does n ot hold f or gradual interfaces 



and domain walls. The opposite limit has been considered in Ref. 19, where Ac = y/2hDo/J (in the paper it was 

£2 

designated by Aj), or with Dq ~ — , Ac ~ £y/2h/ Jr. The limit considered in Ref. implies 2h/ Jt ^ 1 ot Xc ^ i 

T 

and therefore does not hold for ferromagnetic conductors like Fe, Co, Ni and its alloys. 

Semiclassical methods cannot describe processes on length scales smaller than the mean free path, thus cannot 
properly describe abrupt interfaces. It is possible, however, to express boundary conditions in terms of transmission 
and reflection probabilities which ^jconnect the distribution functions on both sides of an interface, and have to be 
computed quantum mechanically.E3 For transport, these boundary conditions translate into interface resistances, 
which arise from discontinuities in the electronic structure and disorder at-the interface. This phenomenon has been 
also extensively studied in the quasi-classical theory of superconductivityEil, where a generalized diffusion approach 
can be used in the bulk of the superconductor, but proper boundary conditions must be used at the interfaces between 
a superconductor and another normal or superconducting metal. j— , 

These effects can be taken into account by flrst-principles band-structure calculations .Ej In collinear systems it is 
possible to circumvent the problem by replacing the interfaces by regions of a flctitious bulk material, the resistances of 
which can be fltted to experiments. This is not possible anymore when the magnetizations are non-coUinear, however, 
because potential steps are essential for the description of the dephasing of the non-collinear spin current and the 
torque. 

We wish to model the multilayer system Fig. || by the diffusion equation and interface boundary conditions. Let 
/(e) be the 2x2 distribution matrix at a given energy e and / the 2x2 current matrix in spin-space. It is convenient 
to expand these matrices into a scalar particle and a vector spin-contribution / — l/o -I- ct • f^, / = (I/q -I- o- • Is)/2. For 
normal metals = If^ -I- ct • s, where is the local charge-related chemical potential and the spin distribution 
function has magnitude /j^ and directions. In the ferromagnet = l/(f + (T-mff = 1(/| +/|)/2 + (t • m(/| — /|)/2, 
where /| and fi are the diagonal elements of the distribution matrix when the spin-quantization axis is parallel to 
the magnetization in the ferromagnet m. 

The diffusion equation describes transport in both the normal metal and the ferromagnet. We flrst consider a single 
interface and disregard spin-flip scattering. The particle and spin currents in the normal metal with diffusion constant 
D are j = Ddxfo — ^^xf^ respectively. The particle and spin currents are conserved: 

i?9|/o^ = 0, D^ff = 0. (1) 
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In the ferromagnet the particle and spin currents are j = D^dxfi + ^I'^xfi and jf = m9a;(D|/| — Difi), where Z)| 
and Di are the diffusion constants for spin-up and spin-down electrons. Current conservation of the spin components 
parallel and antiparallel to the magnetization direction in the ferromagnet read: 

D^dlf^ = Q, Dldlh=0. (2) 

Eqs. are appHcable only inside the bulk layers. The boundary conditions at the interface ariseiom the continuity 
of the particle and spin distribution functions on the normal and the ferromagnetic metal sidesi3'E2l 

fs iN-surface = (/t + /i)/2|F -surface ? (3) 
fj^lN-surface — — /x)/2|F_surface ■ (4) 

Furthermore, particle current is conservedJlS 

[-D9a;/(f ]|N-surface = dx{D^ + -Di/j ) Ip-surface ■ (5) 

We have discussed above why the non-coUinear component of the spin-accumulation decays on a length scale of 
the order of the lattice spacing. This leads to the third boundary condition at the F-N intesfape, namely that the 
spin-current is conserved only for the spin-component parallel to the magnetization direction :eH13 

-surface ■ (6) 



Solving these equatioiiSp-ive recover Eq. ( |A3D with the mixing conductance as found by the magnetoelectric circuit 
theory (Appendix A).Em^l The magneto-electronic circuit theory is thus equivalent to the diffusion approach when 
the system size is larger than the mean free path. However the magneto-electronic circuit theory is a more general 
approach that can also be used for circuits or parts of circuits that are smaller than the mean free path. Note 
that the boundary conditions above do not contain explicit reference to interface conductance parameters and are 
therefore valid only for bulk resistances which are sufficiently larger than the interface resistances. The gain by using 
the diffusion equation is, that we can now easily derive simple analytical results, also in the presence of spin-flip 
relaxation. In the normal as well as ferromagnetic metals spin-flip scattering leads to: 

d 

dxjo = , —is = is/Tsf , (7) 

where the spin-flip relaxation time r,/ is a material dependent parameter. 



III. RESULTS FOR SYSTEMS WITHOUT SPIN-FLIP 



Let us now apply this method to the spin-transfer in a N1-F1-N2-F2-N3 system Fig. | to obtain explicit results for 
the flgure of merit, viz. the ratio of the spin-torque to the charge current through (or voltage bias across) the system. 
The layers are characterized by the lengths Lni, Lfi, Ln2, Lf2, Ln3 and by diffusion constants Dni, Dpi^m), Dn2, 
Dp2,m), Dn3 for each normal and ferromagnetic metal layer, respectively. The resistances of the system are Rni, 
i?Fi,T, Rfi,i, Rn2, Rf2,], Rf2,i and Rn3 with, for example, Rmi = Lni/{AniDmi) and Rpi,] = Lfi/{AfiDfi^]) 
{L and A are length and cross section of a layer respectively). Let us initially disregard spin-flip scattering. 

The continuity of the spin-current at the interface Nl-Fl can easily be shown from Eqs. (0J|,^,^. As a result the 
two layers Nl-Fl behave effectively Hke a single ferromagnetic layer with renormaHzed resistance: 

Rfi,]' = Rfi,] + 2i?Ari , (8a) 

Rfi,i = Rfi,[ + 2Rni . (8b) 

The same is true for the interface F2-N3. As a result it is sufficient to treat only the F1-N-F2 system. In general, 
there are spin-current discontinuities at the interfaces Fl-N and N-F2 which, due to momentum conservation, lead to 
torques acting on the magnetic moments in the ferromagnetic layers. Taking into account all diffusion equations (|]j|) 
and boundary conditions (^,Hj5|J^), and also introducing the parameters R — Rn2, Ri± = {RFi,'\ ± ^-Fj,i)/4, where 
i = 1,2, the torques can be written as 



= iR + R2^){R + R,^)~a^R,^R,^ ^'''^^ - ' (^'^^ 

R + R2+ - aR2-Ri+/Ri- 
^ '°^^-(i? + i?2+)(i? + i?i+)-a^i?i+i?.+ ^"^^ " ' 
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FIG. 2: Torque acting on the first ferromagnet as a function of the relative angle between the two magnetizations for different 
normal metal resistances (the resistances are expressed in units Ri = R2, Pi = 0.5, P2 = 0.2). 



where ti and T2 are torques acting on the magnetizations of the first and second ferromagnet respectively, a 
m2)= cos^, 9 being the angle between the magnetizations. The resistance can also be calculated: 



^{0) = R + Ri+ + i?2+ - 



{1 - a'){Rl_R2+ 



Ri+)/R 



R + R1++ R.2+ + i?i+i?2+(l - 



= (mi 



(10) 



It is worthwhile to rewrite the above formula (^) using the effective polarization P = R-/R+ (which is the polarization 
of a current fiowing through F or N-F layers connected to reservoirs) and the ferromagnet charge current resistance 
Ri = Ri+. The (absolute values of the) torques are then: 



1^2 1 = 



|1 + R/Ri - aPi/P2\ [(1 - a')/oP2] 
(l + i?/i?2)(l + i?/i?i)-a2 

|1 + R/R2 - aP2/Pi\ [(1 - a'VoPi] 

(1 +i?/i?l)(l +i?/i?2) 



(11a) 
(lib) 



As one can see from Eqs. (^) and ( |Tl] ) there is an asymmetry with respect to current inversion. For example, if 
only one polarization can rotate (one ferromagnet is much wider than the other or exchange-biased), domains in 
the two magnetic layers can be aligned antiparallel by currents fiowing in one direction, and reoriented parallel by 
reversing the current fiow. This happens because only one state (parallel or antiparallel) is at equilibrium for a fixed 
direction of the current. If the currents are large enough (depending on other sources of torques such as external fields, 
magneto crystalline anisotropy and damping) the magnetization will fiip, which can be monitored by the change in 
total resistance of 



R{M)-R{]]) 
R{n) 



AR1-R2- 



i?2 + (i?i+ + i?2+)2 - (i?i_ + i?2-) 



(12) 



In the case of unit polarization and R,Ri the relative resistance change j 
predicted by the spin-transfer theorytl and was observed experimentally.!: 



I can be 100 percent. This asymmetry was 
I Note, however, that in these experiments 



the mean free path is comparable to the size of the systems, and the present theory cannot be directly applied. 

From Eq. ( |ll|) follows that the torques are equal to zero for parallel and antiparallel alignments. When the 
numerator of Eq. ( [TTj ) 1 -I- R/Ri[2) — otPi(2) I P2{i) never vanishes, the torque increases with 9 from zero to a maximal 
value which corresponds to an angle larger than 7r/2 and vanishes again when configurations become antiparallel. 
When the nominator of Eq. (||) does vanish for some angle 6*0, the absolute value of torque has a local maximum 
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before (see Fig. g). In principle, it is possible to have an equilibrium magnetization angle 9 = 9q for one current 
direction while equilibrium magnetization angle = or tt for the opposite current direction (this can lead to 
asymmetry for the transition from the anti-aliened state to the aligned state in comparison with the transition from 
aligned to anti-aligned observed experimentallyQ) . 

We propose a setup in which only one magnetization can rotate (usually it is achieved by taking one ferromagnetic 
layer much wider than the other or by exchange biasing) . If one ferromagnetic layer (for example the first one) has a 
resistance Ri <^ R and the other i?2 > R, the torque T2 vanishes whereas the other torque can be simplified to: 

Ti = (l-a2)/oP2 (13) 

The maximal torque in this setup occurs when the magnetizations of the ferromagnet Fl and the ferromagnet F2 are 
perpendicular. 

In general, the spin-torque is maximal when the resistance R of the normal metal vanishes, as could have been 
expected since this also gives the maximum magnetoresistance effect. In Eqs. (^ and ( |ll] ) the size of the magnets 
does not play a dominant role for small normal metal resistances. In this case the torques depend mainly on the 
polarizations. 



IV. RESULTS FOR SYSTEMS WITH SPIN-FLIP 



So far, we have disregarded spin-flip scattering, which can be included readily, though. Here the system N1-F1-N2- 
F2-N3 is analyzed and spin-flip in each normal metal part is considered separately. Introducing spin-flip in Nl and 
N3 leads to the simple result: Eq. (^ without spin-flip remains valid, but with modifled spin resistances: 

sf _ p tanh(Ljvi(jV3)/£,d) 

^ ' \^Nl(NZ)l'^sd) 

where isd is the normal metal spin-flip diffusion length. When L ^ isd, the resistance is governed by the spin-flip 
diffusion length isd, which means that only part of the metal takes part in the spin transfer whereas the rest plays the 
role of the reservoir. This reduction of the active thickness of the device can lead to an effective polarization increase 
by decreasing the effect of Rn in Eq.(^). Spin-flips in the middle normal metal have a larger impact. The torques in 
the presence of spin-flips in N2 read: 

I I n 2^. n P + R'f/Ri-aPi/P2 , 

I I ^^ P + R'^/R2^aP2/Pl . 

|r2| = {1-a KoPi (^^^,,^^^)^^^^,,^^^)_^, , (15b) 
where /3 — cosh{L/£sd) and Pi[2), Ri(2) are given by (^ and (|lj). R^^ is a new effective normal metal resistance: 

R^f = R'jl^^LlM . (16) 

L/£sd 

For L > £sd the torque is significantly reduced by spin-fiips, becoming exponentially small for longer samples. 

Let us now consider spin-flips in the ferromagnet. The treatment of the N1-F1-N2-F2-N3 system is cumbersome, 
so let us concentrate on the simple case of a F-N-F system. In that case formulas remain unchanged, provided i?+ 
and R- are renormalized as: 

„sf _ p tanh(Lfi(f2)/0 ^ 
-^^1(2)- - ^1(2)- 7 77f ' (I'a) 

^Fl{F2)/£sd 



„sf p tanh(Lpi(f2)/-g,d) 

^1(2)+ = ^1(2)+ f 77f • (17b) 

^ ' ^Fl{F2)/£sd 



where £^ is the ferromagnet spin-flip diffusion length. These resistances should be used in Eqs. (^) for the torques in 
F-N-F systems. If spin-flip in the normal metal exists, then the formulas ( p!5| ) should be used. Eqs. ( [1^ ) imply that 
there is no polarization change (as deflned below Eqs. (^)) and only the ferromagnet resistances Ri{2) are affected. 
For L ^ the bulk of the ferromagnet behaves like a reservoir (just like for the normal metal in the same limit) 
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FIG. 3: Torque on each ferromagnet as a function of the relative angle between the two magnetizations for different spin-flip 
diffusion lengths in the normal metal (the resistance R — 0.4 is expressed in units Ri = R2, Pi w P2 = 0.4, L/£sd = 0, 1, 1.5, 
3 the lower plot corresponds to the higher ratio). By dashed line we plotted Slonczewski's resulta for the same polarization. 



and only a slice with thickness £f^ is active. In general, spin-flip in the ferromagnet leads to reduced torques as i?i(2) 
becomes smaller. 

The effect may be quite small as long as the resistance of the ferromagnet is sufficiently larger than that of the 
normal metal (this can also be seen from Eqs. (||) and ([Tl])), so that the polarization of the current is maintained. 

Finally we would like to discuss magnetoresistance and torque for the symmetric Ri+ = R2+ and = i?2-- For 
the angular magnetoresistance we extract from ( p^ the formula observed by Pratt and also shown to be universal for 
any disordered F-N-F perpendicular spin valves in Ref. |2|: 

R{9)-R{0) _ 1- cos 61 

i?(7r)-i?(0) " x(l + cos6') + 2' ^ ' 

with one parameter x that is given by circuit theory 

(19) 

1 — Rer] 

in terms of the normalized mixing conductance 77 = 2g^i/g, the polarization p = igf — 9i)/g, and the average 
conductance g = gi + gi- As we do not take into account interface resistances, in our case the parameters can be 
expressed only via bulk resistances: g = 1/ {R^ + R) + l/{Ri + R), t] = 2/ {Rg), p = 2R- / {2R+ + R). From Eqs. (|l|) 
and ( p^ the analytical expressions of the spin torque on either ferromagnet for current and voltage biased systems 
read: 



' ' x(cos6' + l) + 2 

pg rj sin 9 — fir 

Y(?7-l)cos6'+l + ry 27r 



(21) 



where is the chemical potential in the left (right) ferromagnet. In the presence of spin-flip for the angular 

magnetoresistance we can write (restricting ourself to F-N-F case again): 

m - RjO) ^ (l + x(/3~l)/2)(l-cos0) 

i?(7r)-i?(0) x{P + cos9) + 2 ' ^ ^ 
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where all parameters should be calculated according to ( [l^ ) and ( [l7|) . The dependences of the torque on angle now 
read: 



P(x + l)sing 
(x(cos0 + /3)+2) 



(23) 



P9 



Tysin^ 



where we introduced two parameters: 



2X- 



2 Ai cos^ + ^2 27r 



i 2 Y + 1 



(24) 



(25) 



A. 



X/3- 



X/3 + 2 



Af,sd 



2x + l 



1, K2 



if 



(/?+ !)) + (! + 



1. 



(ki + X^2) 

X + 1 



K2p)(x/5 + 2) 



(26) 



sinh(LAr/^Ar_sd) ' tanh(L_F /£_F srf I 

An interesting result can be drawn from Eqs. ( P4 ) and (23) by comparison with the Eqs. ( |l^ and (po|). In order to 
fit the torque and the magnetoresistance in the presence of spin-fiip we need an additional parameter /? (/3 was defined 
in Eqs. (p!5|)), which depends only on the spin-fiip diffusion length in the normal metal spacer. The general form of 
Eqs. ( p^ ) and ( |2^ ) with only two important parameters seems to be valid even in the presence of interfaces, but this 
has to be confirmed by future studies. Eq. ( p4| ) is cumbersome depending explicitely on the diffusion length in the 
ferromagnets. In Fig. || we have plotted result of Eq. ( [2^ ) for different spin-diffusion lengths in the normal metal. 
The smaller diffusion length corresponds to smaller torques. The curves only qualitatively resemble Slonczewski's 
result for ballistic systems, but it should be pointed out that for the case p = 1, 77 = 2 both approaches result in the 
same formula. 



V. CONCLUSION 



We investigated transport in multilayer systems in the diffusive limit with arbitrary magnetizations in the ferro- 
magnetic layers. The boundary conditions for diffusion equations including spin transfer were discussed and analytic 
expression for the magnetization torques and the angular magnetoresistance were obtained. The torque can be en- 
gineered, not so much via the geometry of the samples (like layer thicknesses) but rather via the materials, the 
ferromagnetic polarization being an important parameter. The asymmetry with respect to the current fiow direction 
has been addressed and the resistance change under magnetization reversal was calculated for different current direc- 
tions. The effect of spin-fiip in the normal metal and ferromagnet was studied analytically. Spin-fiip in the center 
normal metal suppresses the spin-transfer, whereas spin-fiip in the outer normal metals can effectively increase the 
polarization and spin-transfer. The spin-fiip processes in the ferromagnet also diminish the spin-transfer, but not as 
drastically as long as the resistance of the ferromagnet is larger than the normal metal resistance. Finally we show 
in Appendix A that the diffusive approach with carefully chosen boundary conditions leads to results which coincide 
with those from circuit theory. 

We are grateful to Yuli Nazarov, Yaroslav Tserkovnyak and Daniel Huertas-Hernando for stimulating discus- 
sions. This work was supported in part by the NEDO International Joint Research Grant Program "Nano- 
magnetoelectronics", NSF grant DMR 99-81283 and DARPA award No. MDA 972-01-1-0024. 



APPENDIX A: CIRCUIT THEORY APPROACH TO THE DIFFUSIVE SYSTEMS 

Here we show that the diffusion approach is equivalent the circuit r-theory and that the mixing conductance is a 
valid concept also in systems which are dominated by bulk transport.ES We consider a F1-N2-F2 system Fig. |^ (Nl 
and N3 can also be included) connected to two reservoirs Rl and R2 with negligible interface resistances. Note that 
this does not mean that the interface is neglected, because it plays an essential role in the boundary conditions as 
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mentioned in the main text. Since the system is diffusive, a thin slice of a ferromagnet or a normal metal can be 
considered as a node. The mixing conductance can be written in terms of the reflection and transmission coefficients 
and incorporates any kind of details of the contacts, e.g. tunnel, diffusive, and baUistic contacts. We are free to 
define interface resistors via the location of the nodes. Here it is chosen such that the interface width is larger than 
the ferromagnetic decoherence length but smaller than the mean free path. We introduce six nodes: rl is in Rl just 
before the interface, r2 and r3 are before and after the F1-N2 interface, r4 and r5 are before and after the N2-F3 
interface, r6 is in R2 just after the interface Fig. |l|. 

Let us first find the charge and the spin current in Fl (F2) at the interface where the spin-transfer takes place. The 
currents read: 

Jo = (G^ + G^){fi - f^) + (G^ - GT)/f , (Ala) 
I. = m [(G^ - GT)(/o^ - f^) + (G^ + GT)/f ] , (Alb) 

where is the particle distribution function in rl (r6) and Jq is the spin distribution function in r2 (r5). The 
distribution function at r2 and r3 (r4 and r5) is identical due to the continuity boundary condition. The spin-current 
between the ferromagnet reservoir r2 (r5) and the normal metal reservoir r4 (r3) driven by the non-equilibrium 
distributions can be found by using the circuit theory: 

Ii(2) = m2GAr(/f - s • m/f ) - 2ReGTVf (s - (s • m)m) -f (s x m)2ImGTVf , (A2) 

where the spin accumulation in the normal metal reservoir r4 (r3) is given by the unit vector s and the spin distribution 
function . Use was made of = = Gn because r2 (r5) is close to the interface. The component of the 
current perpendicular to the magnetization m is transferred to the magnetization at the interface whereas the parallel 
component is conserved. The torque acting on the magnetization in Fl (F2) therefore becomes: 

Ti(2) = -2ReGl'^/f (s^ - (s^ • m)m) -f (s^ x m)2\mG^^f^ . (A3) 

The mixing conductance is related to the reflection coefficients of an electron from the normal metal to the ferromagnet: 

GM^Y.^5n„.~{rl^rriJ. (A4) 

nm 

Let us now evaluate the mixing conductance for a disordered system. We assume that the junction consists of two 
connected parts. The normal metal section is described by a single scattering matrix for both spin-t and spin- J. 
electrons. The ferromagnetic section requires two independent scattering matrices, one for spin-f and one for spin- J, 
electrons. Scattering at the F-N boundary is disregarded here since it is assumed that the total resistance is dominated 
by the diffuse normal metal and ferromagnetic metal parts of the junction. The total reflection matrix r" for spin-a 
electrons can then be found by concatenating the normal metal and ferromagnetic parts as: 

oo 

r^=rN+ t'^r^F ("-'N^pTtN = rN + x''- (A5) 

71=0 



By inserting ([A5[) into the deflnition for the mixing conductance we flnd that the mixing conductance can be expressed 



as G|| — Gn + SG^i, where 

'^'^Ti = Yl^^^'^ynmXLn] + (xlm (''Jv)"™ )* + ixlm)* xim] ■ (A6) 
nrn 

Eq. ( [A^ ) depends on the phase difference between the scattering paths of spin-up and spin-down electrons. It is 
assumed that there are no correlations between the scattering matrices of the spin-f and spin-| electrons in the 
ferromagnetic part, which is consistent with the a small coherence length. Consequently, in a diffusive systems 
SGii = 0. However, the up- and down-spin parts of the total scattering matrix of the combined normal metal and 
ferromagnetic system are correlated since both spin directions see the same scattering centers in the normal metal 
part. This leads to the conclusion that for diffusive hybrid system: 

Gf, = Gn ■ (A7) 



From Eqs. ( [A1| , [A^JA3| ) and taking into account Eq. ([A7|), and noting that 2GAr = 1/R and Gj = 1/JR|, G| = 1/i? 



one can easily flnd the Eq. (h_ 
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